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Abstract 

Using a single coil mutual inductance technique, we measure the low temperature dependence 
of the magnetic penetration depth in superconducting NbN films prepared with similar critical 
temperatures ~16 K but with different microstructures. Only (100) epitaxial and weakly granular 
(100) textured films display the characteristic exponential dependence of conventional BCS s-wave 
superconductors. More granular (111) textured films exhibit a linear dependence, indicating a 
gapless state in spite of the s-wave gap. This result is quantitatively explained by a model of 
thermal phase fluctuations favored by the granular structure. 

PACS numbers: 74.25.Nf, 74.40. +k, 74.50. +r, 74.80.Bj 



Distinctive property of simple superconducting materials like metals and alloys is the 
absence of thermally activated electronic excitations at low temperature T <C Tc, where Tc 
is the critical temperature. This property is consistent with a model of superconducting gap 
Ak with s-wave symmetry, as predicted by the BCS theory. In contrast, in complex materials 
like heavy fermions cuprates 0] or rhutenates [0], thermal excitations of arbitrarily small 
energies are experimentally observed. This gapless property has been commonly explained by 
an excitation spectrum with nodal points or lines, characteristic of Ak with p- or d-symmetry 
respectively. According to theoretical considerations, non-s symmetries are favored by non- 
phonon mediated pairing mechanisms like antiferromagnetic spin fluctuations @] . However, 
experimental and theoretical studies show that a gapless state can be induced in s-wave 
superconductors by residual interactions coexisting with the pairing interaction, such as 
magnetic impurities 0], spin density waves 0, proximity effects H or phase fluctuations 
of the order parameter |^ |10], 0. In conclusion, a non-s symmetry implies a gapless state, 
but the opposite does not hold. 

In this work we investigate whether the above predictions of gapless state apply to the 
s-wave superconductor NbN. By measuring the variation of the magnetic penetration depth 
A(T) — A(0) = SX{T) at T > 1.5 K, we study the low energy excitations on films with 
different microstructures obtained by varying growth conditions. NbN is a model system 
for our purpose taking into account the following: i) the simple cubic structure with fee 
lattice |T^|; ii) the conventional BCS s-wave superconducting properties [T^l; Hi) the rela- 
tively high critical temperature Tc ~ 16 K enables to study the low energy excitations in 
a wider temperature region; iv) various materials preparation routes have been developed 
and optimized, especially in view of electronics applications |l^, |15|, |16 . 



For this work we have selected six NbN films prepared by DC magnetron sputtering with 
similar close to the bulk value f«16 K but with different degrees of granularity obtained 
by varying growth conditions. These conditions are reported in Table I (see ||15| , [l^ for more 
details): growth temperatures are either room temperature or 620 °C; substrates are (1102) 
AI2O3 or (100) MgO single crystals or (1102) AI2O3 single crystals with a 12 nm thick (100) 
MgO buffer layer; NbN thicknesses range from 10 to 1400 nm. Three different types of 
microstructure are reproducibly obtained: i) (100) epitaxial; ii) highly (100) textured; Hi) 
(111) textured. Using (100) MgO buffer layers, a high degree of (100) texture is obtained 



only in sufficiently thin {d < 150 nm) NbN films |]T6 
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The film microstructure is characterized by X-ray diffraction in the Bragg- Brent ano ge- 
ometry, X-ray pole figure measurements and transmission electron and atomic force mi- 
croscopies (TEM, AFM). X-ray rocking curve analysis indicates that samples A1043 and 
A1013 are epitaxial with FWHM of the (420) Bragg refelection less than 1°. The former 
sample is thin (10 nm), whereas the latter is thick (950 nm). The film on the MgO buffer 
layer (A1044) is thin (10 nm) and highly (100) textured. The remaining three films on 
sapphire (A1057, A1060 and A1063) are all thick (280-1400 nm) and display predominant 
(111) texture. The different degrees of granularity in these films are confirmed by AFM and 
TEM. For example, in sample A1043, we find an atomically flat surface |l^ and a crystal 
structure with apparently no grain boundaries (see fig. la), which indicates that this film is 
monocrystalline. A granular structure is evident in the (111) textured sample A1060 from 
the TEM cross-sectional micrograph of fig. lb. We note a columnar structure with grains 
of size 15-35 nm and well-defined grain boundaries, as previously reported for films grown 



under the same conditions [|T5 



The microstructural properties of the six films are correlated with DC resistivity g and 
DC transport critical current density jc measurements. The results of these measurements 
are summarized in Table I. In all films, T^. is close to the bulk value ~16 K, which indicates 
that the microstructure does not appreciably perturb the superconducting phase. Smaller 
residual resistivities f)(20 K) 50/xf2cm and larger resistivity ratios f)(300 K)/f)(20 K) ^ 1.2 
are found in the epitaxial films. The highest values of jc (4.2 K) > 3 MA/cm^ are found 
in the epitaxial or thin (100) textured samples, while the three thick (111) textured ones 
exhibit values < 1 MA/cm^ that decrease with increasing NbN thickness. In agreement with 
earlier studies [|l^ and with our TEM results, we conclude the following: i) the epitaxial 
samples have higher crystalline quality; ii) in the (111) textured samples, the degree of 
granularity - defined as the degree of grain decoupling - increases with increasing thickness. 

To correlate the above structural and transport properties with the low energy excitations, 
we have carried out 5A(T) measurements in the range 1.5 K-Tc using a single-coil mutual 



inductance technique described elsewhere |T8[. Typical frequency and magnitude of the e.m. 
field perpendicular to the film surface are 2-4 MHz and < 0.1 mT respectively. Care has been 
taken to always operate in the linear response regime, i.e. below the Josephson critical field 
-Bel, monitored by varying the applied magnetic field. Note that we measure an effective 
magnetic penetration depth that, in the dirty limit, differs from the London penetration 
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depth by a correction factor of the order of unity [0. This hmit is appropriate for our 
films, for which we estimate a mean free path ^1 nm and a coherence length ^ r:^ 3 — 5 nm 
T5| , p!7| . As a good approximation, the above correction factor is temperature independent in 



the limited temperature range of our measurements. Therefore the temperature dependence 
of the variations SX{T) is not affected by dirty limit effects. 

In figs. 2a-b, we report the experimental SX{T) curves obtained. Only the epitaxial and 
thin (100) textured films display the exponential behavior 6X ~ exp(— A/A;bT), character- 
istic of BCS 5-wave superconductors (see fig. 2a). All of the three thick (111) textured 
samples exhibit a linear behavior which can not be approximated by any power law (see fig. 
2b); the range of linearity is limited below ^ 2.7 K in the thinner film A1057 and extends 
up to ~4 K in the thicker films A1060 and A1063. A BCS s-wave fit perfectly accounts for 
the experimental data of the three samples of fig. 2a. A(0) is the only free parameter of this 
fit, whereas the ratio 2A(0)/A;b^c is kept fixed at 4.6, taken from tunneling data at 1.5 K 
on sample A1043. The fitted A(0) values are in the range 200-400 nm, in agreement with 



previous reports [T^, |T^]. A BCS fit is satisfactory also for the thinner (111) textured film 
A1057 of fig. 2b; indeed, the behavior of this sample is intermediate between exponential 
and linear. For the remaining two samples A1060 and A1063, no BCS fit is possible for any 
plausible set of parameters. 

These results give evidence of a gapless superconducting state in the (111) textured films. 
We recall the mechanisms that may induce a gapless state in a conventional BCS s-wave 
superconductor: i) magnetic impurities; ii) proximity effects; in) phase fluctuations of the 
order parameter. In our case, the flrst two possibilities are excluded: i) the purity of the Nb 
target used is 99.995 % and traces of magnetic impurities are below 10 ppm. We conclude 
that the amount of these impurities in the fllms is negligible. In any case, the gapless samples 
contain no more impurities than the others, since all of them are prepared with the same Nb 
target; ii) a previous TEM and tunneling study indicates that the intergranular phases 
are insulating, thus no proximity effects are expected. The third picture remains plausible: 
as suggested earlier [^, |10|, [T^, ^ |22|, both quantum and thermal phase fluctuations can 
be relevant in granular superconductors with short coherence length ^ like our NbN fllms. 

To verify the validity of the above picture in our case, we consider an array of super- 
conducting grains with intergranular Josephson coupling energy Ej and charging energy U, 



as flrst proposed by Abeles EM, where U is associated with the excess of Cooper pairs n 
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in each grain. U and Ej are related to the intergranular capacitance C and the critical 
current Ic respectively, according to the relations U = 2e^/C and Ej = hlc/2e. Following 
previous work we assume the array to form a square lattice of constant a equal to 

the average grain size. The resulting Hamiltonian is: 

H=^Y.^^ + EjJ2[l- cos(^. - ^,)] (1) 

i i^j 

where is the phase of the order parameter in the i-th grain and rii = —2id/ddi is its 
conjugate operator. 

Eq. (1) describes both quantum and classical phase fluctuations. Two arguments suggest 
a classical description of our results in the range T > 1.5 K of our experiment, i) According 
to theoretical studies [y, |23| supported by numerical calculations , quantum fluctuations 



lead to a reduction of the renormalized superfiuid density, implying a progressive decrease of 
dX/dT as T ^ K. This effect is not observed in the curves of fig. 2b which remain linear 
above 1.5 K. ii) In our (111) textured films, we estimate U ^ Ej, which corresponds to a 



purely classical regime 0, |2^, Indeed, typical figures are: a ^ 20 nm, A(0) ~ 300 nm, 
jc ~ 0.5-1 MA/cm^, thickness and relative dielectric constant of the intergranular regions 
respectively t < 2 nm and > 2. Thus, we obtain Ej ^ 50-100 meV and U < 10 /leV. 

In the classical regime, the fluctuation strength is determined by the ratio Ej/ksTc, where 
Tc is the mean fleld critical temperature. Using the above estimate, we flnd Ej ^ ksT^ in 
our granular fllms. Thus, according to previous studies p2| , ^5] , we expect small fluctua- 
tions which do not inhibit long range phase order at any temperature. Indeed, we do not 
observe any appreciable decrease of Tc induced by granularity (Tc is even slightly larger in 
the (111) textured fllms). These fluctuations coexist with the usual quasiparticle excitations 
predicted for BCS s-wave superconductors. Both excitations lead to a progressive increase 
of A as temperature increases. In the Gaussian approximation, valid at low temperatures 
T <C Ej/ks and in the absence of dissipation, only long wavelength longitudinal fluctu- 
ations of arbitrarily small energy are created (Goldstone modes). This produces a linear 
increase 5A ~ T [P, |10|, [ill experimentally observed in flg. 2b. On the other hand, the cre- 
ation of quasiparticles requires a flnite thermal energy > A, which leads to the exponential 
dependence of flg. 2a. 

We now analyze in quantitative form the experimental curves of flg. 2b within the frame- 
work of the above picture. Using Table I, we estimate flrst the fluctuation contribution and 



then the temperature at which the quasiparticle contribution becomes significant. Following 
previous calculations in the Gaussian regime, our lattice model with U=0 leads to: 



dX 
df 



ekB 1 



(2) 



where Ic ~ jc^A at K and we neglect the weak temperature dependence of this product 
at low temperature. In Table I we report the values of dX/dT estimated using eq. (2) 
and compare them with the experimental values. Taking into account the uncertainty in 
the estimate of a, we notice a good agreement. Interestingly, the trend of jc in the (111) 
textured films follows the trend of dX/dT, as predicted by eq. (2): the thicker the film, the 
smaller jc and the larger the slope dX/dT. In the films of fig. 2a, exhibiting a BCS s-wave 
behavior, is larger. Thus, the fiuctuation effect vanishes and only the quasiparticle effect 
is observed. To estimate the temperature at which the two effects become comparable, we 
separately evaluate their contribution to dX/dT using Table I and the BCS fit of fig. 2b. 
We find that the linear increase of A produced by thermal phase fiuctuations is expected 
to dominate at T < 4 K. This is in agreement with the point of departure from linearity 
in samples A1060 and A1063. As discussed before, in sample A1057 a gapless behavior is 
visible but does not dominate yet above 1.5 K, since a BCS fit is still satisfactory. 

In conclusion, we have reported for the first time experimental evidence of gapless state 
in a conventional BCS s-wave superconductor, NbN. Evidence is found in granular films 
exhibiting a clear linear dependence of A in the range 1.5-4 K. Epitaxial or less granular 
films exhibit the conventional dependence of s-wave superconductors and no trace of gapless 
state is found. A simple model of granularity- induced thermal phase fiuctuations of the order 
parameter (Goldstone modes) quantitatively accounts for the above linear dependence. No 
signature of quantum fiuctuations is observed above 1.5 K, as expected taking into account 
the relatively large grains of our films. 

Our evidence of granularity- induced gapless superconductivity in NbN raises the question 
whether this phenomenon is observable in other granular superconductors. In cuprates, 
which display gapless properties in agreement with a d-wave symmetry model, it should be 
investigated if these properties could be - at least partially - explained by intrinsic granular 
properties associated with chemical or electronic inhomogeneities in the coherence length 
scale < 3 nm. In any case, our results suggest that the only rehable methods for testing 
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non-s symmetries are those based on phase-sensitive measurements. 
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FIG. 1: a). Cross-sectional TEM image of the (100) epitaxial sample A1043. Note the monocrys- 
talline structure evidenced by the contrast produced by the diagonal (110) planes. The notation 
"R" indicates the (1102) orientation of the sapphire substrate, b). The same as in a) for the thick 
(111) textured A1060 sample. Note the granular structure made of columnar grains. 

FIG. 2: a). Temperature dependence of the magnetic penetration depth A measured in the films 
with conventional BGS s-wave behavior: A1013 (full circles), A1043 (squares) and A1044 (open 
circles), b). The same as in a) for the films with gapless behavior: A1057 (lozenges), A1060 (full 
triangles) and A1063 (open triangles). The BCS-like curve of sample A1044 (circles) of a) and the 
corresponding BGS s-wave fit (solid line) are also shown for comparison. 



TABLE I: Summary of sample properties and comparison between experimental values of dX/dT 
and theoretical ones obtained using eq. (2). For the samples following the BGS thermally activated 
behavior, A(0) is obtained using a BGS s-wave fit, as discussed in the text. 



sample 


A1043 


A1013 


A1044 


A1057 


A1060 


A1063 


NbN thickness [nm] 


10 


950 


10 


280 


680 


1400 


substrate 


AI2O3 


MgO 


AlaOg/MgO 


AI2O3 


AI2O3 


AI2O3 


substrate orient. 


(1102) 


(100) 


(1102)/(100) 


(1102) 


(1102) 


(1102) 


growth temp. [°G] 


620 


room temp. 


620 


room temp. 


room temp. 


room temp. 


microstructure 


epitaxial 


epitaxial 


highly textured 


textured 


textured 


textured 


NbN orientation 


(100) 


(100) 


(100) 


(111) 


(111) 


(111) 


Tc[K] 


14.2 ±0.1 


15.5 ±0.09 


14.4 ±0.1 


13.4 ±0.1 


16.1 ±0.1 


16.3 ±0.1 


jc [MA/cm^] 


3.0 lb 0.7 




9.0 ±0.3 


1.3 ±0.7 


1.1 ±0.3 


0.64 ± 0.04 


q{20 K) [nQcm] 


53 


50 


147 


448 


88 


67 


q{300 K)/g{20 K) 


1.2 


1.2 


0.9 


0.7 


0.9 


1.0 


A(0) [nm] (BGS fit) 


240 


260 


270 


470 






5A/5Texp. [A/K] 











2.1 ±0.5 


3.3 ±0.2 


5.6 ±0.5 


dX/dT th. [A/K] 











1.6 


2.2 


3.3 
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